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5_ G MSK SIGNAL PROCESSORS F OR TMPROVED rOMMUNICATTONS 

CAPACITY AND QUALITY 



BACKGROUND OF THE INVENTION 

1. Field of the Invention 

10 The present invention pertains generally to reducing interference in wireless 

communications systems using Gaussian Minimum Shift Keyed (OMSK) or other 
Minimum Shift Keyed (MSK) signals, and more specifically to a method and 
apparatus for (i) rejecting interfering signals received at a communications 
receiver and occupying the same or adjacent spectral band as one or more desired 

15 signals, (ii) separating, if needed, multiple desired communication signals received 
at a communications receiver, and (iii) correcting distortion that affects one or 
more desired communication signals received at a communications receiver such 
as that caused by a multipath radio-frequency (RF) propagation channel. 

2. Description of the Background Art 

20 The background of the present invention falls into two general categories: 

(i) work pertaining to the filtering structures that perform interference rejection, 
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signal separation, and distortion correction; and (ii) work pertaining to the 
adaptive algorithms that adjust the coefficients in thes'e filtering structures to allow 
them to operate in signal environments that are endemic to mobile, portable, and 
personal wireless communication systems. This work, however, differs 
substantially from the present invention. 
A. Filtering Structures 

With the exception of the single-sensor FRESH filter and the Joint 
Maximum Likelihood Sequence Estimator, all prior work on interference rejection 
and signal separation that could be applied to a system that receives multiple 
GMSK or MSK signals in the same spectral band an<^tefier interval requires that 
multiple sensors (antennas) be used. Another requirement of this body of work is 
that the spatial characteristics of any two signals received by the multiple sensors 
must be sufficiently different; if this requirement is not met (or is poorly met) then 
the filtering structure will not be able to reject the interference and separate the 
desired signals (or will do so with poor quality). This body of work includes all 
filtering structures that belong to the class of linear time-invariant (LTI) spatio- 
temporal filters; this class of filters is defined by the fact that their outputs can be 
modeled mathematically as sums of amplitude-scaled, phase-shifted, and time- 
shifted versions of the output signals of the multiple sensors. Notably, this class, 
of filters characterizes virtually all filtering structures that have been previously 
developed to address the problem (and problems similar to it) solved by the 
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present invention. 

With regard to the single-sensor FRESH filter, none of the prior works (i) 
discuss or describe application of FRESH filtering to GMSK or other MSK signal 
types, (ii) discuss or describe application to GSM or DECT standards, (iii) suggest 
using frequency shifts equal to plus or minus one-half the bit rate of a signal or 
equal to two times the carrier frequency plus or minus one-half the bit rate, (iv) 
show a 3-path LCL-FRESH structure relevant to multiple GMSK signals with 
equal carrier frequencies and bit rates which (for complex envelope with zero 
carrier offset) contains no unconjugated frequency-shifted paths, (v) show a 2- 
path, 2-input LCL-FRESH structure relevant to adjacent-channel GMSK 
interference, (vi) show a 3-path, 3 -input LCL-FRESH structure relevant to 
adjacent-channel GMSK interference, (vii) show a multi-sensor 3-path LCL- 
FRESH structure comprising one each of identical versions of a 2-path or 3-path 
LCL-FRESH structure at the output of each sensor with all^3-path LCL-FP^SH 
filter outputs summed, (viii) describe any partially blind methods for adaptation of 
any filter to separate GMSK or MSK signals, (ix) describe any training-aided 
constant modulus algorithm (CMA) methods for adaptation for any signal type, 
(x) demonstrate separability of any GMSK or MSK signals using any FRESH 
filtering method, (xi) describe the theoretical spectral redundancy properties of 
GMSK that provide for signal separability using FRESH filtering, or (xii) show 
the measured spectral redundancy properties of GMSK signals. Furthermore, the 
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single-sensor FRESH filter has a severe limitation that is transcended by the 
present invention: it cannot accommodate multiple sensors in any way that allows 
it to reject and/or separate more than 2 co-channel and 2 adjacent-channel GMSK 
or MSK signals, in contrast to the present invention. 

Although a Joint Maximum Likelihood Sequence Estimator (JMLSE) 
operating with a single sensor may theoretically have the ability to reject and/or 
separate more than 2 GMSK or MSK signals, the quality of the separated signals 
must degrade substantially when the number of these exceeds two. Furthermore, 
the computational complexity of the JMLSE is extremely high and is not practical 
to implement in real-time without expensive hardware that draws copious amounts 
of electrical energy. These drawbacks become more severe (at an exponential 
rate) as the number of sensors to be accommodated increases beyond one and the 
number of signals to be jointly separated increases. 

The filtering structures of the present invention do not suffer from the 
severe implementational drawbacks of the JMLSE; they accommodate single or 
multiple sensors to great advantage; and they can reject and/or separate multiple 
signals even when one or more pairs of these signals cannot be rejected and/or 
separated by the LTI spatio-temporal filters. More specifically, the filtering 
structures of the present invention provide capabilities, for rejecting interfering 
signals and separating desired signals and correcting their distortion, that are not 
provided by any prior filtering structures. 
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B. Adaptation Methods 

The filtering structures used in this invention, and in other means and 
apparatus for interference rejection, signal separation, and correction of distortion 
are controlled through the selection of values of filter parameters or coefficients. 
In some embodiments of these filtering structures, these coefficients can include 
the gains and phases used to weight various signals prior to summing. In typical 
wireless communication systems, the values of these coefficients must vary 
properly because the spatio-temporal characteristics of the received signals vary in 
time. That is, the values of the coefficients of these filtering structures must be 
chosen adaptively. 

With the exception of the present invention, there are two classes of 
methods for adaptively choosing these coefficient values: (i) conventional 
methods that rely exclusively on the knowledge of training sequences embedded . 
in the desired signals, this knowledge being used to choose an initial setting for 
the coefficient values, typically according to a least-squared error criterion; and 
(ii) blind methods that do not use embedded training sequences but instead exploit 
one or more properties of the interference or of the desired signals. 

Conventional methods that use a known training signal are well-known to 
those skilled in the art of signal processing, and are described in the open literature 
pertaining to adaptive signal processing and adaptive filtering. Conventional 
methods are limited in their applicability, however, because the number of 
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coefficient values that they can adapt successfully must be less (much less for 
good performance) than the number of known samples within the embedded 
training sequence. This number of known samples is fixed by the wireless 
communication system (e.g., for all signals in GSM that convey voice traffic, this 
5 number is the product of 26 and the number of effectively independent samples 
per bit (which is typically 2) used in the digitization of the continuous-time 
received signals, assuming a discrete-time implementation). With regard to the 
problem addressed by the present invention, in which the filtering structures can 
have, potentially, a large number of adjustable coefficients, the number of known 

10 samples within the embedded training sequence can be insufficient to provide a 
useful initial setting. A filtering structure initialized in this way would not provide 
sufficiently high quality estimates of the desired signals to be useful in the 
communications system. 

In contrast, blind methods can make use of the entire record of the received 

15 signals over which. the spatio-temporal characteristics (i.e., angles of arrival, 
relative delays, and relative phases of the multiple propagation paths of the 
^ multiple signalsjand carrier and bit phases of the transmitted signals) are 

^ approximately constant. This circumvents the limit^conventional methods on 

the number of filter coefficients that can be adapted. However, the various blind 

20 methods have other severe drawbacks that prevent them from operating 
successfully in many situations, including the one addressed by the present 
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invention. Specifically, the class of Constant Modulus Algorithms (CMA) and 
related methods can require much longer data records, to converge to a high- 
quality solution, than are available. Furthermore, their convergence to a solution 
is not guaranteed, leading to a confusion between interfering signals (which 
should be rejected) and desired signals (which should be separated from each 
other and from the interference). Another class of blind methods relies on 
estimating the directions of arrival of the received signals; these estimates are then 
used in a table lookup (through the array manifold, which is a table of calibration 
data for the multiple-sensor receiver) to find the coefficients of an LTI spatio- 
temporal filter. However, these so-called direction- finding methods tend to 
perform poorly in the presence of severe multipath, unknown number of signals, 
and unknown interference and noise fields, these three impairments being endemic 
to the wireless communication systems of interest in the present invention. 
Furthermore, the direction- finding methods do not provide the types of parameter 
estimates that are needed to adapt the filtering structures of the present invention. 

The adaptation methods of the present invention exploit in a novel manner 
both the embedded training sequences and the constant modulus property of the 
desired signals, and in doing so circumvent the undesirable aspects of the 
conventional methods and of the blind methods. That is, the adaptafion methods 
of the present invention, which choose the values of coefficients in filtering 
structures, such as those of the present invention, provide capabilities that are not 
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provided by prior adaptation methods. 



SUMMARY OF THE INVENTION 



The present invention provides a method and apparatus fo r rejectin g 
undesired interference and separating and corr ecting the distortion of desired 

5 interfering si gnals of the Gaussian Minimum Shift Keyed (GMSK) or other 
Minimum Shift Keyed (MSK) type, which arise for example in wireless 
communication systems based on either the Digital European Cordless 
Telecommunications (DECT) standard or the Global System for Mobile 
Communications (GSM) standard and, as a result, improving the quality and 

10 quantity of service offered by a wireless communications system that uses GMSK 
or other MSK signals. 

The utility of the present invention may be appreciated, for example but not 
exclusively, by manufacturers of wireless communications equipment and 
providers of wireless communications services who wish to improve the quality of 

15 communications offered to users of their systems, and/or who wish to increase the 
number of users of their systems (e.g., to increase revenue in a commercially 
operated system). 

Quantity of service improves when a base station is able to support a greater 
number of users for a fixed frequency allocation by, for example, being able to 
20 reduce frequency reuse distance, extend range, possibly reuse frequencies within a 
cell, and more frequently assign users to adjacent channels. Quality of service 
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improves, for example, when requests made to a central coordinator (hereafter 
referred to as the base station) by users to initiate calls are rejected less frequently 
(i.e., blocking probability is reduced), when calls in progress by mobile users are 
not abnormally terminated (i.e., probability of dropped calls is reduced), and when 
the intelligibility of the voice messages in the call is better (i.e., the parties to the 
call hear less noise, interference, abnormal silences, and the like). 

The present invention addresses these general problems by means of a 
method and apparatus that perform the following functions: (1) rejection of 
interfering signals received at a communications receiver and occupying the same 
(or adjacent) spectral band as one or more desired signals, (2) separation (if 
needed) of multiple desired communication signals received at a communications 
receiver, and (3) correction of distortion that affects one or more desired 
communication signals received at a communications receiver (e.g., distortion 
caused by a multipath radio-frequency (RF) propagation channel). 

Accordingly, an object of the invention is to provide an array of one or more 
antennas whose outputs are processed by conventional RF front-ends and are then 
processed by a filtering structure that rejects interference, separates one or more 
desired signals, and corrects distortion to obtain high-quality estimates of desired 
signals. 

Another object of the invention is to provide a class of filtering structures 
that can reject interfering signals, separate one or more desired signals, and correct 



9 



EM519050807US 




the distortion in these desired signals. 

Another object of the invention is to provide a class of filtering structures 
that can operate on the outputs of any number of antennas after these outputs are 
downconverted by conventional RF front-ends, including the special case when 
5 only one antenna is used. 

Another object of the invention is to provide a linear-conjugate-linear 
(LCL) frequency-shift (FRESH) filter that uses the particular pattern of frequency- 
shifts and conjugations to exploit the cyclostationarity properties that are 
characteristic of GMSK and other MSK signals, thereby providing the desired 
10 interference-rejection, signal-separation, and distortion-correction capabilities. 

Another aspect of the invention is to provide a linear-conjugate-linear 
polyperiodically time-variant (LCL-PTV) filter that uses the particular pattem of 
conjugations and time- variant filter coefficients to exploit the cyclostationarity 
properties that are characteristic of GMSK and other MSK signals, thereby 
15 providing the desired interference-rejection, signal-separation, and distortion- 
correction capabilities. 

Another object of the invention is to provide an LCL-FRESH fractionally- 
spaced equalizer (LCL-FRESH-FSE) that operates on higher-than-bit-rate- 
sampled data to provide bit-rate-sampled output and uses the particular pattem of 
20 frequency-shifts and conjugations to exploit the cyclostationarity properties that 
are characteristic of GMSK and other MSK signals, thereby providing the desired 
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interference-rejection, signal-separation, and distortion-correction capabilities. 

Another object of the invention is to provide an LCL-FSE filter that 
operates on higher-than-bit-rate-sampled and frequency-shifted (by one-fourth of 
the bit rate) data and uses the sum of FSE processed data and conjugated data to 
exploit the cyclostationarity properties that are, characteristic of GMSK and other 
MSK signals, thereby providing the desired interference-rejection, signal- 
separation, and distortion-correction capabilities. 

Another object of the invention is to provide a real-imaginary (RI-FSE) 
filter that provides a computationally efficient alternate embodiment of the LCL- 
FSE filter. 

Another object of the invention is to provide an LCL filter that operates on 
bit-rate-sampled and frequency-shifted (by one-fourth of the bit rate) data and uses 
the sum of linearly filtered data and conjugated data to exploit cyclostationarity 
properties that are characteristic of GMSK and other MSK signals, thereby 
providing the desired interference-rejection, signal-separation, and distortion- 
correction capabilities. 

Another object of the invention is to provide a real-imaginary (RI) filter that 
provides a computationally efficient alternate embodiment of the LCL filter. 

Another object of the invention is to provide a class of methods for adapting 
the adjustable coefficients in any of the filtering structures of this invention so 
that the filter can properly perform the fiinctions of rejecting interference, 
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separating one or more desired signals, and correcting distortion. 

Another aspect of the invention is to provide a means of adapting the 
adjustable coefficients, by using both the known training sequence(s) embedded in 
the desired signal(s) and the knowledge that the desired signal(s) have constant 
modulus, thus providing training-augmented constant-modulus (TACM) 
adaptation. 

Another object of the invention is to provide a means of adapting the 
adjustable coefficients, by using the known training sequence(s) embedded in the 
desired signal(s) to constrain adaptation that uses the knowledge that the desired 
signal(s) have constant modulus, thus providing training-constrained constant- 
modulus (TCCM) adaptation. 

These general aspects and their details give this invention significant 
capabilities of great utility in wireless communication systems. The capability 
provided by this invention can be used to improve the quality and quantity of 
communication services provided by such systems. The quality of service can be 
improved by using the invention to suppress interference in a particular frequency 
band and time slot (jointly referred to as a channel) and spatial cell caused by 
other users in the same or adjacent channels and in the same or adjacent spatial 
cells and to correct signal distortionjjand the quantity of service can be increased 
by reallocating channels in a manner that intentionally results in more interfering, 
but separable by this invention, signals per channel and cell, in order to 
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accommodate more users. 

These improvements are made possible by this invention because this 
invention, using only one sensor, can separate multiple spectrally overlapping 
signals of certain types and can correct even severe frequency-selective fading and 
other distortion and, because the number of such signals that can be separatedby 
this invention, using multiple sensors increases proportionately with the number of 
sensors. In cellular communication systems that use base stations, the allocation 
of channels can (but need not) be made so that the introduced interference occurs 
only in up-link channels (from mobile units to base stations) so that the invention 
need only be implemented in the base station receivers and not in the mobile 
receivers. The invention can also be used at base stations (or at mobile units, or 
both) using only a single receiving antenna or multiple antennas. 

With reference to the specific example of GMSK and other MSK signals, 
between 2 and 4 GMSK and other MSK signals in a single channel can be 
accommodated with one antenna, depending on how many signals arise from users 
in the same channel and how many arise from adjacent channels. More generally, 
with a number Mof antennas, between 2Mand 4M signals can be accommodated. 

Further objects and advantages of the invention will be brought out in the 
following portions of the specification, wherein the detailed description is for the 
purpose of frilly disclosing preferred embodiments of the invention without 
placing limitations thereon. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be more fully understood by reference to the following 
drawings which are for illustrative purposes only: 

FIG. 1 is a graph showing the cross-coherence functions of a sample GMSK 
5 signal and its conjugated and frequency-shifted versions. 

FIG. 2A is a graph showing the spectral content of a GMSK (0.3) signal, 

x(t). 

FIG. 2B is a graph showing the spectral content of a conjugated and 

frequency shifted version of the signal shown in FIG. 2A, y^{t) = x* (t)e~-^'^-^^^ . 
10 FIG. 2C is a graph showing the spectral density of the degree to which x 

shown in FIG. 2 A andj^y shown in FIG. 2B are correlated. 

FIG. 3A is a graph showing the spectral content of a GMSK (0.3) signal, 

x(t). 

FIG. 3B is a graph showing the spectral content of a conjugated and 

15 frequency shifted version of the signal shown in FIG. 3 A, yjit) = x*{t)e^'^^^^ . 

FIG. 3C is a graph showing the spectral density of the degree to which x 
shown in FIG. 3 A and j;^ shown in FIG. 3B are correlated. 

FIG. 4 is a functional block diagram of a 3 -path LCL-FRESH filter in 
accordance with the present invention for performing interference rejection and/or 
20 co-channel separation and/or distortion correction for GMSK or other MSK 
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signals. 

FIG. 5 is a functional block diagram of an LCL-PTV filter in accordance 
with the present invention for performing interference rejection and/or co-channel 
separation and/or distortion correction for GMSK or other MSK signals. 

FIG. 6 is a functional block diagram of a "type- 1" LCL-FPIESH-FSE filter 
in accordance with the present invention for performing interference rejection 
and/or co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 7 is a functional block diagram of a "type-2" LCL-FRESH-FSE filter 
in accordance with the present invention for performing interference rejection 
and/or co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 8 is a functional block diagram of a "type- 1" LCL-FSE filter in 
accordance with the present invention for performing interference rejection and/or 
co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 9 is a functional block diagram of a "type-2" LCL-FSE filter in 
accordance with the present invention for performing interference rejection and/or 
co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 10 is a functional block diagram of a "type-1" lU-FSE filter in 
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accordance with the present invention for performing interference rejection and/or 
co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 1 1 is a functional block diagram of a "type-2" RI-FSE filter in 
accordance with the present invention for performing interference rejection and/or 
co-channel separation and/or distortion correction for GMSK or other MSK 
signals. 

FIG. 12 is a functional block diagram of an LCL filter in accordance with 
the present invention for performing interference rejection and/or co-channel 
separation and/or distortion correction for GMSK or other MSK signals. 

FIG. 13 is a functional block diagram of an RI filter in accordance with the 
present invention for performing interference rejection and/or co-channel 
separation and/or distortion correction for GMSK or other MSK signals. 

FIG. 14 is a functional block diagram of a 5-path LCL-FRESH filter in 
accordance with the present invention for performing interference rejection and/or 
co-channel^eparation and/or distortion correction, with improved cancellation of 
adjacent channel interference^for GMSK or other MSK signals. 

FIG. 15 is a functional block diagram of a 2-input, 2-path LCL-FRESH 
filter in accordance with the present invention for canceling adjacent upper 
channel interference. 

FIG. 16 is a functional block diagram of a 2-input, 2-path LCL-FRESH 
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filter in accordance with the present invention for canceUng adjacent lower 
channel interference. 

FIG. 17 is a functional block diagram of a 3-input, 3-path LCL-FRESH 
filter in accordance with the present invention for canceling both adjacent upper 
and^ower channel interference. 

FIG. 18 is a functional block diagram of a joint space-time multi-sensor 
FRESH filter in accordance with the present invention. 

FIG. 19 is a functional block diagram of space-time factorized processor in 
accordance with the present invention. 

FIG. 20 is a data structure timing diagram for a time slot in a GSM system, 
showing the number and location of each type of bit. 

FIG. 21 is a data structure timing diagram showing the time periods Tq and 
for the time slot shown in FIG. 20. 

FIG. 22 is a functional block diagram of a receiver apparatus in accordance 
with the present invention. 

FIG. 23 is a functional block diagram of an embodiment of the coherent M- 
antenna collection apparatus block shown in FIG. 22. 

FIG. 24 is a functional block diagram of an embodiment of the filter 
apparatus block shown in FIG. 22. 

FIG. 25 is a functional block diagram of an embodiment of the frequency- 
shifting, time-shifting, and polarized switching apparatus block of the filter shown 
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in FIG. 24 where the filter is configured as a 3 -path LCL-FRESH filter. 

FIG. 26 is a functional block diagram of an embodiment of the frequency- 
shifting, time-shifting, and polarized switching apparatus block of the filter shown 
in FIG. 24 where the filter is configured as a type-1 LCL-FRESH-FSE filter. 

FIG. 27 is an altemative embodiment of the frequency-shifting, time- 
shifting, and polarized switching apparatus shown in FIG. 26 where the filter is 
configured as a type-2 LCL-FRESH-FSE filter. 

FIG. 28 is a fimctional block diagram of an embodiment of the frequency- 
shifting, time-shifting, and polarized switching apparatus block of the filter shown 
in FIG. 24 where the filter is configured as a type-1 LCL-FSE filter. 

FIG. 29 is an altemative embodiment of the frequency-shifting, time- 
shifting, and polarized switching apparatus shown in FIG. 28 where the filter is 
configured as a type-2 LCL-FSE filter. 

FIG. 30 is a fimctional block diagram of an embodiment of the frequency- 
shifting, time-shifting, and polarized switching apparatus block of the filter shown 
in FIG. 24 where the filter is configured as a type-1 RI-FSE filter. 

FIG. 3 1 is a fimctional block diagram of an embodiment of the frequency- 
shifting, time-shifting, and polarized switching apparatus block of the filter shown 
in FIG. 24 where the filter is configured as an LCL filter. 

FIG. 32 is a fimctional block diagram of an embodiment of the frequency- 
shifting, fime-shifting, and polarized switching apparatus block of the filter shown 
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in FIG. 24 where the filter is configured as an RI filter. 

FIG. 33 is a functional block diagram depicting the polarized switch portion 
of the frequency-shifting, time-shifting, and polarized switching apparatus shown 
in-FiGr^through FIG. 32. 

FIG. 34 is a functional block diagram of an embodiment of the tapped time 
delay (TDL) apparatus block of the frequency-shifting, time-shifting, and 
polarized switching apparatus shown in FIG. 25 through FIG. 32, where the input 
and output vectors can be real or complex. 

FIG. 35 is a functional block diagram of a TACM adaptation apparatus in 
accordance with the present invention. 

FIG. 36 is a functional block diagram of a TCCM adaptation apparatus in 
accordance with the present invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring more specifically to the drawings, for illustrative purposes the 
present invention is generally shown with reference to FIG. 1 through FIG. 36 
where like reference numerals denote like parts. It will be appreciated that the 
apparatus may vary as to configuration and as to details of the parts and that the 
method may vary as to the steps and their sequence without departing from the 
basic concepts as disclosed herein. 

A, Filtering Structures and Their Theoretical Foundations 

The filtering structures of the present invention exploit the 100% spectral 
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redundancy exhibited by GMSK and other MSK signals. That is, these filtering 
structures exploit the fact that information conveyed in one subband of the GMSK 
or other MSK signal is also conveyed in another subband. In fact, for any GMSK 
or other MSK signal, the waveforms in any two such related subbands are nearly 
perfectly correlated. Thus, this spectral redundancy augments the more widely 
known and exploited spatial redundancy in which the waveforms due to a given 
source and received by two different antennas in an array are perfectly correlated. 

Just as a receiver without multiple antennas and the spatial filtering 
structure needed to properly combine their outputs loses out on the well-known 
performance gains due to the exploitation of spatial redundancy, a receiver 
without the filter structures of the present invention loses out on the significant 
performance gains due to the exploitation of spectral redundancy. 

By exploiting both spatial redundancy and spectral redundancy (also 
referred to as spectral correlation), the filtering structures of the present invention 
provide at least twice the level of processing gain of an LTI spatio-temporal 
filtering structure. For example, an M-antenna array can separate up to 2Mto 4M 
GMSK or other MSK signals using these filtering structures, or can separate up to 
the more conventionally attainable number Mof GMSK or other MSK signals but 
can do so with significantly higher reliability and quality than a conventional 
spatial filter. These capabilities can be explained by reference to the spectral 
correlation of the GMSK or other MSK signals, and by reference to the 
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relationship between correlation (regardless of its origin or type) and linear 
filtering. 

B. Relationship Between Correlation and Linear Filtering 
Although the relationship between correlation and linear filtering is 

5 fundamental, it is not widely understood even by those skilled in the art of signal 
processing. Therefore, for purposes of describing the present invention we will 
explain this relationship by a simple, although abstract, analogy: if 
x(t) = a(t) + b(t), and y(t) = c(t) + d(t), where a(t) and c(t) are highly correlated, 
and where b(t) and d(t) are highly correlated, then an appropriate linear 

10 combination of x(t) and y(t) (i.e., the output of an appropriately chosen linear 
filtering structure applied jointly to x(t) and y(t)) can provide high-quality 
estimates of one or more of the signals a(t), b(t), c(t), and d(t). 

C. Spectral Correlation of GMSK and Other MSK Signals 

To understand how such performance can be possible, and to understand 
15 what filtering structures are applicable to GMSK and other MSK signals, it is 
necessary to understand the spectral correlation properties of GMSK and other 
MSK signals. For this purpose, these properties are best described by the cross- 
coherence function, which is a frequency-domain cross-correlation coefficient. 
Specifically, let x(t) denote the complex envelope of a GMSK signal having zero 
20 carrier- frequency offset and bit rate and \Qiy(t) be a conjugated and frequency- 
shifted version of x(t), where the frequency-shift is either fjl or -/y2: 
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. y,{t) = x\t)e-j^''(f>'^'>' (1) 
y^{t) = x\t)e'^^-(f>'^^)' (2) 

W\\hy(t) =y,(t) OYy(t) =y2(t). We consider the cross-coherence function C^y(0 
defined by 

CM)= , (3) 

where Sjfy(f) denotes the cross spectral density between x(t) and y(t), and and 
Syy(f) denote the power spectral densities ofx(t) and y(t), respectively. The 
magnitude \C^(f)\ versus frequency /for a GMSK signal is shown in FIG. 1, where 

the solid line represents C^^ (/) and the dashed line represents C^y^ (/) , and is 



clearly very nearly equal to one over the entire frequency channel (-100 kHz to 
+100 kHz) allocated to the signal. 

For any / such that \C^y(f)\ is nearly equal to one, the narrowband 
components of x(t) and y(t) at frequency / are nearly perfectly correlated. For our 
definitions ofx(t) and y(t), this is equivalent to stating that the narrowband 
component of x(t) at / and the narrowband component of x*(t) at either / + f^2 or 
/ - fi/2 are nearly perfectly correlated. Referring to FIG. 2A through FIG. 2C, the 

fact that (/) in FIG. 2C is nearly equal to unity for narrowband components 
numbered -6,-5, ...,-1,0 implies that each pair of like-numbered narrowband 
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components taken from x(t) in FIG. 2 A and y^{t) = x*(/)e"^^*^*^in FIG. 2B are 
nearly perfectly correlated. The same is true for the continuum of pairs of 
narrowband components in x(t) andyjft) at identical values of / Accordingly, an 
appropriate linear combination of one of these components from x(t) and the 
corresponding one from yj(t), can provide the same level of processing gain 
against other signals, interference, and noise as an appropriate linear combination 
of the signals received by at least twice as many antennas as used to receive x(t). 
Similar statements apply to linear combinations involving x(t) and 

J^zCO = X* {t)e^^^^^ . Thus, a filtering structure that uses linear combinations of 
x(t),yj(t), and y2(t) can provide the same level of processing gain against other 
signals, interference, and noise as an appropriate linear combination of the signals 
received by twice as many antennas as used to receive x(t). 

Based on these facts about the spectral correlation properties of GMSK and 
other MSK signals, two specific interference scenarios are cited here as examples 
to illustrate, but not to limit, the applicability of the filtering structures of this 
invention to practical problems in wireless communication systems. For 
simplicity, only a single antenna is used in both scenarios; the capability of the 
proposed filtering structures increases proportionately with the number of 
antennas used. 
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Co-Channel Interference Rejection 
A single antenna receives a desired GMSK signal and a co-channel GMSK 
signal having the same carrier frequency and comparable signal power. Neither 
signal can be recovered by a linear time-invariant filtering structure. However, 
with reference to FIG. 2A through FIG. 2C, each of the components numbered 
-6 to 0 (and the components in the continuum from -150 kHz to 0 kHz) in x(t) can 
be linearly combined with its corresponding component in y^ft) to yield a clean 
estimate of the desired signal while rejecting the co-channel interference. 
Similarly, with reference to FIG. 3 A through FIG. 3C, each of the components 
numbered 0 to 6 (and the components in the continuum from 0 to 150 kHz) in x(t) 
can be linearly combined with its corresponding component in y2(t) to yield a 
clean estimate of the desired signal while rejecting the co-channel interference. 
Thus, by using both yj(t) and y 2(0, the desired signal can be recovered over the 
entire band from -150 kHz to +150 kHz, or over just the allocated GSM frequency 
channel from -100 kHz to +100 kHz if desired. In this way, the filtering structures 
in this invention can reject a GMSK co-channel interfering signal even when only 
one antenna is used. 

Ad jacent Channel Interference Rejection 
A single antenna receives a desired GMSK signal and an adjacent channel 
GMSK signal possibly having much greater signal power. Specifically, let the 
adjacent channel signal have a carrier frequency that is 200 kHz (per GSM carrier 
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allocation protocols) less than that of the desired signal. Due to spectral leakage, 
the adjacent channel interferer degrades the desired signal, which cannot be 
recovered by a linear time-invariant filtering structure. However, with reference 
to FIG. 3 A through FIG. 3C, each of the components numbered 4 to 8 (and the 
components in the continuum from 100 kHz to 200 kHz) in the adjacent channel 
interference (now downconverted to zero carrier offset), which comprise the 
spectral leakage into the desired signal, can be linearly combined with its 
corresponding component in yj(t) to cancel out the effects of the spectral leakage. 
Similarly, an adjacent channel interferer whose carrier is 200 kHz above that of 
the desired signal can be canceled out. In this way, the filtering structures in this 
invention can reject two GMSK adjacent channel interfering signals even when 
only one antenna is used. 

Signal Distortion Correction 
When a transmitted signal undergoes propagation through a non-ideal 
channel, such as a multipath propagation channel, the received signal can be a 
severely distorted version of the transmitted signal, including for example 
frequency-selective fading (i.e., deep spectral nulls). However, with reference to 
FIG. 2 A through FIG. 2C, any of the narrowband components numbered -6, -5, . . . 
-1,0 (and the components in the continuum from -150 kHz to 0 kHz) in x(t) that 
have been severely attenuated by frequency-selective fading (to such an extent 
that they cannot simply be amplified without strongly increasing the noise level), 
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can be replaced with corresponding properly phase-adjusted and amplitude-scaled 
components 'myi(t), thereby removing the spectral nulls. Similarly, with reference 
to FIG. 3 A through FIG. 3C, spectral nulls in the band from 0 kHz to +150 kHz ' 
can be removed. 

D. Specific Filtering Structures 

Various filtering structures in accordance with the present invention will 
now be described with reference to FIG. 4 through FIG. 17. In the filtering 
structures depicted therein, all blocks with only a letter inside denote standard 
linear time- invariant filters (ignoring time variation due only to adaptation), those 
blocks with a letter followed by (t) denote ^nearly periodically time-variant filters, 
and those blocks with a letter and FSE denote linear time-invariant fractionally 
spaced equalizers (filters), all of which are known s ymbols - used in the signal 
processing art. Furthermore, blocks with the symbol ( )* inside denote 
conjugators, blocks with "Re" inside denote real-part selectors (extractors), and 
blocks with "Im" inside denote imaginary-part selectors (extractors), each of 
which operates on a complex signal and each of which is also referred to herein as 
a polarized switcher. These elements are known to those having ordinary skill in 
the signal processing art or are as otherwise described herein. 

L CL-FRESH Filter 

The preceding observations about \C^y(f)\ and the role of high or nearly 
perfect correlation in motivating linear filtering structures are taken into account 
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in the 3-path time-domain linear-conjugate-linear (LCL) frequency-shift filter 
(FRESH filter) 1 0 in accordance with the present invention whose structure for 
GMSK and other MSK signals is depicted in block-diagram form in FIG. 4. When 
applied to the signal x(t) received at a single antenna, two copies of this structure 
or the other filter structures described below can separate and/or correct the 
distortion of two co-channel GMSK or other MSK signals. The 3-path LCL- 
FRESH filter is appropriate for data that is sampled at a rate highei; than the b it 
rate. Accordingly, the outp ut of th e 3-path LCL-FRESH filter must be sampled at 
the bit rate prior to being processed, for example, by a GMSK demodulator that 
operates on bit-rate sampled data to recover the bit stream. 

LCL-PTY Filter 

As an altemative to the 3-path LCL-FRESH filter, the effect of the two 

complementary frequency shifts e-^^^^^ and e"*^^*^*^ and the associated filters g_ and 
g+ can be achieved also with a single time-variant filter g(t). The resulting 
structure is referred to herein as a linear-conjugate-linear periodically time- variant 
(LCL-PTV) filter 12 shown in block-diagram form in FIG. 5. The filter g(t) has a 
time-variant impulse response according to the possibly periodic time-variations 
in the possibly cyclostationary statistics of the received signals. In some 
applications, it is appropriate for the coefficients of time-variant filter g(t) to vary 
with frequency yj,, in which case there exists a one-to-one mapping between the 
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LCL-PTV filter and the 3 -path LCL-FRESH filter. In other applications, such as 
those in which the time-slot boundaries of the signals do not line up sufficiently 
well and the length of a time-slot is a non-integer number of bit periods (e.g., as in 
GSM), the coefficients of filter g(t) must change every time one or more of the 
received signals encounters a time-slot boundary. The LCL-PTV filter is 
appropriate for data that is sampled at a rate higher than the bit rate. Accordingly, 
the output of the LCL-PTV filter must be sampled at the bit rate prior to being 
processed, for example, by a GMSK demodulator that operates on bit-rate sampled 
data to recover the bit stream. 

LCL-FRESH-FSE Filter 
In both the 3-path LCL-FRESH filter and the LCL-PTV filter, the sampling 
rate of the filter output exceeds the bit-rate (typically for a linear filter processing 
GMSK, it is an integer multiple of the bit-rate yj,, such as in which case the 
filter output is typically sampled at the bit-rate prior to being processed by the 
demodulator (which provides an estimate of the bit stream conveyed in the desired 
signal). Thus, an altemative to the 3-path LCL-FRESH filter is to move this bit- 
rate sampling operation back through the filtering structure to the extent possible 
without changing the capability of the filter and to eliminate redundant operations. 
Altemative embodiments of the resulting filter structure 14 (type-1) and 16 (type- 
2), which are referred to herein as LCL-FRESH fractionally-spaced equalizers 
(LCL-FRESH-FSE), are depicted in block-diagram form in FIG. 6 and FIG. 7, 
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respectively. This filter differs significantly from the LCL time-invariant filter by 
its particular use of conjugated signal paths and by its special use of frequency 
shifts to properly exploit the cyclostationarity properties of the GMSK or other 
MSK signals. The bit-rate-sampled signal estimate s{k) must be processed by a 



filter to obtain a filter structure in which the output of each bit-rate-sampler is 
frequency-shifted by This shifted signal is then split into two paths: the upper 

10 path is processed by an LTI filter, and the lower path is conjugated and processed 
by another LTI filter. That is, the lower path following each bit-rate-sampler no 
longer contains a frequency shift of (-1)*". However, the output final summer must 
be frequency-shifted by /. Altemative embodiments of the resulting filter 
structure 18 (type-1) and 20 (type-2) are shown in FIG. 8 and FIG. 9, respectively. 

15 This filter structure will be referred to herein as an LCL-FSE filter rather than an 
LCL-FRESH-FSE because no path is frequency-shifted any differently than any 
other path. 

The LCL-FSE filter can be simplified without any loss of performance by 
observing that the two paths at the output of each bit-rate sampler and 
20 frequency-shifter form an LCL filter, which can be equivalently realized by 
modifying the upper path to operate on the real part (instead of operating on the 
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demodulator to recover the bit stream. 



L CL-FSE and RI-F$E Filters 



Straightforward manipulations can be performed on the LCL-FRESH-FSE 
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entire complex signal) and by modifying the lower path to operate on the 
imaginary part (instead of operating on the conjugated complex signal). Since this 
structure replaces the two LCL filters by two Real-Imaginary (RI) filter pairs, the 
resulting structure is referred to as the RI-FSE filter and can be interpreted to be an 
alternate but mathematically equivalent embodiment of the LCL-FSE filter. 
Altemative embodiments of this filter structure, 22 (type-1) and 24 (type-2), are 
shown in FIG. 10 and FIG. 11, respectively. 

LCL and RI Filters 

The LCL and RI filters are related to the LCL-FSE and la-FSE filters, but 
operate on bit-rate-sampled data. This is possible for GMSK, and for some other 
MSK signals, because the bit-rate is larger than the RF bandwidth occupied by the 
signal, and so the bit-rate satisfies Nyquist's sampling criterion for this signal. 

The result of simplifying the LCL-FSE filter as a result of the input being 
sampled at the bit-rate is the LCL filter 26, shown in FIG. 12. 

The LCL filter can be simplified without any loss of performance by 
observing that it can be equivalently realized by modifying the upper path to 
operate on the real part (instead of operating on the entire complex signal) and by 
modifying the lower path to operate on the imaginary part (instead of operating on 
the conjugated complex signal). Since this structure replaces the LCL filter by a 
Real-Imaginary (RI) filter pair, the resulting structure 28, shown in FIG. 13, is 
referred to as an RI filter and can be interpreted to be an alternate but 
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mathematically equivalent embodiment of the LCL filter. 

The RI filter offers the same implementation advantage over the LCL filter 
as the RI-FSE offers over the LCL-FSE: since the data paths being processed by 
various LTI or FSE filters are real, both the application of the filters to the data 
5 and the adaptation of the filters require less arithmetic than if the data paths were 
complex. Significantly, since the RI filter is mathematically equivalent to the 
LCL filter, there is no loss of quality in the estimated signal at the output. 



LCL, and RI filters are described herein as if their filter coefficients are fixed over 
one time-slot, they can in fact be made time-variant (as with the LCL-PTV filter) 
so as to accommodate the changes in statistical properties of the received signals 
when, for example, one or more of the signals encounters a time-slot boundary. 



^ 15 However, in systems in which there is E^tgn ^ degrcc o f control over the 

synchronization of the time-slot structures of the multiple co-channel signals, it is 
advantageous to force the time-slot boundaries of the multiple signals to be as 
close together as possible. Such desired synchronization could certainly arise if 
the multiple signals originate from mobile units operating in the same cell in a 
20 wireless cellular communication system; similarly, if the multiple signals originate 
fi-om mobile units operating in different cells, then the desired synchronization 



Time-Variant Versions of LCL-FRESR LCL-FRESH-FSE. 
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Although the 3-path LCL-FRESH, LCL-FRESH-FSE, LCL-FSE, RI-FSE, 
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could be made to arise through cooperation of the base stations in these different 
cells. 

Variations for Improved Rejection of Adjacent Channel Interference 
The aforementioned filtering structures are ideally suited for the rejection of 
co-channel interference and separation of co-channel desired signals and removal 
of signal distortion. However, variants of these filtering structures can provide 
rejection of adjacent channel interference by exploiting the spectral correlation 
properties of the adjacent channel interference in addition to those of the co- 
channel signals. Capability to reject adjacent channel interference can be 
especially important when the adjacent channel signals are much stronger than the 
desired channel signals (e.g., due to the near-far problem). The aforementioned 
filtering structures can be generalized (or modified) to exploit the cyclostationarity 
properties of the adjacent channel interference in addition to those of the desired 
channel signals (or, if desired, only those of the adjacent channel interference). 
For example, the 3-path LCL-FRESH filtering structure can be generalized to 
include two additional paths, one for each of the upper and lower adjacent 
channels that can be present. The signals present in these two additional paths are 
highly correlated with the interference, from the corresponding adjacent channels, 
that contaminates the desired signals; thus, correctly linearly combining these 
additional paths with the other three paths can significantly reduce the adjacent 
channel interference. The resulting 5-path LCL-FRESH filtering structure 30 is 
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depicted in FIG. 14, wherein xj(t) and x^(t) are the complex envelope signals from 
the lower and upper adjacent channels (i.e., downconverted to complex baseband 
from their respective RF carriers), respectively. 

LCL-FRESH filtering structure 30 shown in FIG. 14 can be simplified as sh0wn in 
FIG. 15 through FIG. 17 by deleting the two paths feeding into filters g/and g_ 
can be deleted, yielding a new filtering structure, with only three pafths, whose 
functionality is limited to the cancellation of adjacent channeUnterference, instead 
of joint cancellation of adjacent- and co-channel interfereilce, signal separation, 
and distortion correction. These two-input two-path ^Ix^lL-FRESH filtering 
structures 30a, 30b, which are shown in FIG. /fsMl-LCL-FRESH-U) and FIG. 16 
(2 I-LCL-FRESH-L), respectively, can be u^d to remove interference from either 
the upper adjacent channel by using inpdt x^(t) and frequency shift If^ - iJ2^ or the 
lower adjacent channel by using imxit x_(t) and frequency shift + 4/2, 
respectively. As an altemativ^ interference from both upper and lower adjacent 
channels can be simultangfously removed using the 3-input 3-path LCL-FRESH 
filter 30c shown in pen ? (3 I -LCL- F R ES II-UL'). 

Based on the aforementioned descriptions of the manner in which the LCL- 
PTV, LCL-FRESH-FSE, LCL-FSE, RI-FSE, LCL, and RI filtering structures can 
be obtained from the 3-path FRESH filter, those skilled in the art of signal 
processing will appreciate that these filtering structures can also be generalized to 
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improve their ability to reject ad jacent channel interference. Such generalization 
may be accomplished by altering these filtering structures to exploit the 
cyclostationarity properties of the adjacent channel interference. 

Finally, it is noted that combinations of the outputs of various filtering 
structures may be used advantageously. For example, each of two copies of the 3- 
path LCL-FRESH filter can be applied to each of the lower and upper adjacent 
channels to provide improved estimates of the adjacent channel signals; these " 
estimates can then be appropriately frequency-shifted, adaptively filtered, and then 
subtracted from the received data in the desired frequency channel, thus reducing 
the adjacent channel interference. The signal resulting thereby can then be " 
processed by the 3-path LCL-FRESH filter to separate co-channel signals and 
correct for distortion. The adjacent channel interference canceling filters can 
either be adapted jointly with the co-channel separation and distortion correction 
filters or be adapted separately. 

Filter Generalization 

It will be appreciated that each of the foregoing filter structures include 
filters (comprising time-shifters and linear combiners), fi-equency shifters, 
polarized switchers (conjugators or real- or imaginary-part selectors), and 
summers (linear combiners). Thus, common core elements of all of these filtering 
structures are time-shifters, frequency shifters, polarized switchers, and linear 
combiners. Table 1 below maps conjugate cycle frequencies and input and output 
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sampling rates into frequency-shift values and polarized switch functions (e.g., 
conjugation or real- or imaginary-part extraction) defining the specific filter 
structures described herein. 

-^Fable-t- 



Conjugate Cycle Frequencies, Input/Output 
Sampling Rates 


Filtering Structures ^y^^^^ 


±4/2, 0/0 


FIG. 4, LCL-FRmi 
(FIG. 5, LCLtFTV) 


±fh/2, 0/B 


FIG. 7^CL-FRESH-FSE-2 
(n0:6, 8, 9, 10, 11) 


±4/2, B/B / 


FIG. 12, LCL 
(FIG. 13) 


2f, - fb/2, 0/0 


FIG. 15,2I-LCL-FRESH-U 


-2f, + V2, 0/0 


FIG. 16, 2I-LCL-FRESH-L 


2f, - fb/2, -2f, + iyjiy^ 


FIG. 17, 3I-LCL-FRESH-UL 


±fi/2, -li^yXli, - fb/2, 0/0 


FIG. 14, 5-pathLCL-FRESH 



5 / 

whey^g^ = the carrier frequency of the lower adjacent channel, f^.' = the carrier 
i^equency of the lower adjacent channel, O = oversampled, and B - bit rat e 
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• sampled . The filt e r structu r es listed i n parenthesis^;ejnathemaH O{rfly-equi yatent 
to the basic structure no^^^^ntKesIvand the former can be derived from the 
l^te£. ^rfitandard -bkTck diagram manipulations. 

It will further be appreciated that each of the foregoing filtering structures 
are baseband structures designed for complex baseband signals corresponding to a 
desired down-converted channel centered at zero fi^equency, and if the desired 
channel is centered at a non-zero fi"equency, then standard baseband-to-passband 
transformations can be applied to these baseband filtering structures to obtain their 
passband counterparts for processing passband signals. 

Joint Space-Time Filtering and Factorize d Filtering 

As depicted in FIG. 18 and FIG. 19, any of the aforementioned filtering 
structures can be used in either of two ways: (1) in a joint space-time filtering 
structure 32, in which each antenna output is processed by one of these filters, and 
then the filtered outputs are summed, or (2) in a space-time factorized filtering 
structure 34, in which a purely spatial filter is applied to the received array data, 
and then the spatially filtered data is processed by one of these filters. The former 
structure has the capability to attain higher performance at the expense of greater 
complexity than the latter. Alternatively, the 3-path LCL-FRESH filter can be 
replaced by either the LCL-PTV filter or the 5-path LCL-FRESH filter in either 
FIG. 18 or FIG. 19, or by any of the LCL-FRESH-FSE, LCL-FSE, RI-FSE, LCL, 
and RI filters filter in either FIG. 18 or FIG. 19 to yield a bit-rate-sampled estimate 
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s{k) of the desired signal waveform. 

Variations: Single vs. Multiple Sensors, 
& Single vs. Multiple Desired Signals 

By using the spectral redundancy theory of cylcostationary signals, we have 
proven mathematically that the best performing method for combining FRESH- 
filtering capability and spatial-filtering (multi-sensor linear combining) capability 
is to replicate the best single-sensor FRESH filter by passing the data from each 
sensor through its own FRESH-filter prior to linearly combining the multiple data 
sets. This results in a special form of periodically time-variant linear-conjugate 
spatio-temporal filter. For example, using the single-sensor 3 -path LCL-FRESH 
filter with M sensors yields the 3-path M-sensor LCL-FRESH (3-path MS-LCL- 
FRESH) spatio-temporal filter. Therefore, it follows that the LCL-FRESH, LCL- 
PTV, LCL-FRESH-FSE, LCL-FSE, RI-FSE, LCL, and RI filtering structures can 
apply not only to single-sensor receivers (in which case x(t) is a scalar signal and 
the LTI or FSE filters in each signal path are single-input filters) but also to 
multiple-sensor receivers (in which case x(t) is a vector signal and the LTI or FSE 
filters in each signal path are multiple-input filters). Similarly, the LCL-FRESH, 
LCL-PTV, LCL-FRESH-FSE, LCL-FSE, RI-FSE, LCL, and RI LCL filtering 
structures can be applied to estimate a single desired signal (in which case the LTI 
or FSE filters in each signal path are single-output filters; in this case, if multiple 
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signals are desired to be estimated, then multiple copies of the filtering structure 
can be implemented in parallel) as well as to estimate multiple desired signals (in 
which case the LTI or FSE filters in each signal path are multiple-output filters). 
It will also be appreciated that the frequency-shifters and conjugators employed in 
the filters, as well as the real and imaginary part extractors, become multiple-input 
devices. These observations emphasize the fact that the novelty and utility of the 
filtering structures reside largely not in the dimensions of the signals in any 
particular signal path within the structure but rather within the patterns of 
conjugations, real-part and imaginary-part selections, frequency shifts, and 
sampling operations and the ways in which linear combinations of signals 
subjected to these transformations are formed. 

Implementation of LTI and FSE Filters 
It will be appreciated by those skilled in the art of signal processing that the 
linear-time-invariant (LTI) filters and FSE filters describing the various filtering 
structures of the present invention can be implemented in numerous ways. In one 
such implementation, an LTI filter is implemented as an FIR filter^and is applied 
to its input by using the well-known overlap-and-add or overlap-and-save 
algorithms that utilize the fast Fourier transform (FFT). In another such 
implementation, an LTI filter is implemented using a tapped delay line, multipliers 
attached to each tap, and a summer. In yet other implementations, the LTI filter 
might be implemented using an infinite-impulse response (IIR) filter instead of an 
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FIR filter; however, it is noted that adapting such filters can be computationally 
expensive and/or unreliable. 



Finally, the interaction between the filtering structures (used for interference 



5 rejection and/or co-channel signal separation) and the demodulator (used for 
recovering the bit stream from the estimate of the desired signal) must be 
considered. Several demodulators are currently in use for GMSK signals: the 
maximum likelihood sequence estimator (MLSE) based on the Viterbi algorithm, 
the multiple-symbol differential detector, and the decision feedback equalizer 

10 (DFE). Any one of these demodulators can be used following the filtering 

structures to demodulate the estimated signal of interest. Altematively, since the 
DFE is itself a filtering structure, it is reasonable to merge the demodulation 
function implemented by the conventional DFE and the interference-rejection 
and/or signal-separation functions implemented by the space-time filtering 

15 structure; the resulting filtering structure would have the desirable property of 
very low complexity relative to the receivers based on MLSE or multiple-symbol 
differential detection, and may be capable of attaining a higher-quality recovered 
bit stream due to its smaller number of adjustable coefficients for the same 
memory length. 

20 Because all of the altemative filtering apparatus described herein include 

the capability of the feedforward part of the DFE, the result of merging the 
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Ct^ filtering apparatus with the DFE is that ono aimply follows the filtering apparatus, 

•^l; - be it single- or multiple-sensor, witly the feedback part of the DFE. Thus, the 
feedback part alone functions as the demodulator. Algorithms for adaptively 
adjusting this demodulator are not within the scope of this invention. 
5 Similarly, the complexity of the MLSE demodulator can be reduced and its 

performance can possibly be increased by the effects of the distortion-removal 
capability of the FRESH filtering structures. For example, the requirements of the 
channel estimator that is an integral part of the MLSE demodulator can be eased 
by virtue of the distortion removal or even just distortion reduction. 
10 D. Adaptation Methods 

The following discussion of adaptation methods is made in the context of 
the GSM system, which is one type of wireless communication system to which 
the present invention applies. Those skilled in the art of communication systems 
and signal processing will appreciate that these adaptation methods can apply also 
15 to any other communication system that uses signals with embedded training 
sequences and constant modulus. 

In GSM, each time slot lasts 577 |is (156.25 bit periods at 270.833 
kbits/sec) and the bits therein are structured as shown in FIG. 20. The three tail 
bits on either side of the main burst are all zeros, and the guard time lasting 8.25 
20 bit periods prevents consecutive time slots from colliding due to different bulk 
propagation delays. The 1 16-bit message in each slot can be either digitized 
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speech (after vocoding, convolutional coding, and interleaving) or control 
information. The 26-bit training sequence embedded in the middle of the slot is 
used in existing conventional GSM receivers, for example, to perform channel 
estimation prior to MLSE demodulation. 

In the present invention, not only is this convenient training sequence 
exploited but also the very usefiil constant-envelope (or constant-modulus) 
property of GMSK and other MSK signals is exploited to adapt the filter 
coefficients. Four adaptation algorithms are explained, the first two of which are 
conventional and the second two of which are novel and unique to the present 
invention: 

1 . conventionally trained adaptive equalization; 

2. constant-modulus (CM) blind adaptive equalization; 

3. training-augmented constant-modulus (T ACM) partially blind 
adaptive equalization; and 

4. training-constrained constant-modulus (TCCM) partially blind 
adaptive equalization. 

The first two adaptation algorithms provide the building blocks for the second 
two. All four of these algorithms are directed to finding the values of the filter 
coefficients in a matrix W that controls a linear combiner (the number of columns 
in this matrix is set to be equal to the number of desired signals to be estimated 
and the number of rows is set to be equal to the number of weights in the linear 
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combiner to be adaptively adjusted). The linear combiner applies these complex 
weights to the elements of a VQCtov y(t) containing various differently delayed, 
frequency-shifted, and possibly conjugated versions of the received data from the 
M-sensors (denoted by the vector x(t)) to produce a vector of signal estimates 

5 s{t) = W y(/) , where superscript H denotes conjugate transposition. For 
example, in the 3 -path LCL-FRESH filter, let each of the three filters be 
implemented by a tapped delay line (TDL) containing L taps, followed by a linear 
combiner. Then the signal estimate s(t) can be expressed as 



s(t) = 



x(t) 
x(t - 1) 



x{t-L^ 1) 
X {t)e ^ 



X it-\)e 



X {t-L + l)e 0 
X {t)e ^ 



X {t-X)e 



X {t-L^-X)e 



(4) 



10 where y(t) is defined in the obvious way. Similar definitions of y(t) can be 



EM519050807US 



constructed for the LCL-FRESH-FSE, LCL-FSE, RI-FSE, LCL, and RI filtering 
structures, as well as for the 5 -path LCL-FRESH filter and any extensions of 
these in which the number of paths, patterns of conjugations or real-part and 
imaginary-part selectors, choices of frequency shifts, and origins of input signals 

5 for the various paths are altered. Thus, in the following descriptions, no further 
reference to the specific dependencies ofy(t) onx(t) is made, since all of the 
filtering structures of the present invention can be represented using this common 
framework. Furthermore, to simplify the descriptions only the case of W having a 
single column is considered; to emphasize that the matrix degenerates to a vector 

10 in this case, the set of filter coefficients is denoted by w. 

Conventionally Trained Adaptive Equalization 
In conventionally trained adaptive equalization, the receiver knows that a 
training signal So(t) is transmitted by the desired user starting at time toi and ending 
at time to2. This knowledge and the received data are used to find w. Specifically, 

15 w is chosen so as to minimize the time-averaged squared error between the 
training signal So(t) and the estimate 5(/)over the training interval (i.e., for all 




(5) 



where O^^^ denotes time-averaging over the period tQ^<t<tQ2, and the 
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correlation matrix between any two vectors (or scalars) a(t) and b(t) is defined by 



time interval containing the embedded training signal. This method is simple to 
implement, with low computational complexity when implemented using the 
5 recursive least squares (RLS) algorithm or any of the various fast algorithms (e.g., 
the modular multichannel method or the fast affme projection algorithm), and it 
converges reliably and quickly when the length of the training sequence exceeds 
the number of filter coefficients (i.e., the length of w). However, in the presence 
of strong co-channel interference and/or when the length of w is comparable to or 



^ 10 greater than the^rarning'signal, this conventionally trained adaptive equalization 



algorithm may not provide sufficiently reliable or high quality signal estimates. 
Constant Modulus (CM^ Blind Adaptive Equalization 
In total contrast to the conventionally trained adaptive equalizer, the 
constant-modulus (CM) blind adaptive equalizer ignores completely any training 
15 signal information. Instead, in this application it exploits the fact that OMSK and 
other MSK signals have constant modulus (constant envelope), a property which 
is degraded or destroyed by the presence of co-channel interference and channel 
distortion. The CM algorithm attempts to solve the following minimization 



Rib ^ = (a(Ob^(0) • With reference to GSM systems and FIG. 21, is the 




problem: 
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min( 



^ k(Ol 



(6) 

^ where^^^n our case is the set of all time instants within a particular time slot. 

Since this is a nonlinear minimization problem with no closed-form solution, and 
a least-mean-squares (LMS) stochastic gradient descent algorithm would converge 
5 much too slowly, an iterative block approach (the so-called least-squares CM 
algorithm) is used. At the J^^ iteration, the weight vector is found by solving 

^^"^ min(|wfy(0-^*-i(0f) (7) 



^ H 

where z^_i (0 = ^^-'^^ , . That is, the modulus normaHzed term in the squared 



wf-iy(0|j 

error is obtained from the filter output during the previous iteration. This allows a 
10 simple linear solution to be found: 

These two steps, finding the modulus-normalized filter output and solving this 
line'ar equation, are iterated until convergence. The CM algorithm allows the 
equalizer to be adapted successfully in the absence of a training signal, although it 
15 is known to exhibit unreliable convergence in some situations and its reliability 
may be lower than desired when the observation interval is short (as it is for GSM 
time slots). Another complication with the CM algorithm arises when it is applied 
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to environments containing multiple co-channel signals. Multi-target CM 
algorithms have been developed to allow multiple sets of filter weights to be 
adapted, each to extract a different signal. A subsequent signal sorting operation 
must then be performed to identify the correspondence between extracted signals 
and users; this correspondence can be achieved easily if the users are assigned 
different training sequences. Notice though that this training signal information 
might be more profitably used for adaptation of the filter weights themselves 
rather than simply to sort the output of the blindly adapted filters. This 
observation is pursued through two other adaptive algorithms, discussed below. 

Training- Augmented C M (TACM) 
Partially Blind Adaptiv e Equalization 
Unlike the conventionally trained adaptive equalizer and the CM blind 
adaptive equalizer, each of which ignores significantly usefiil information that is 
exploited by the other, the TACM partially blind adaptive equalizer exploits both 
the training signal available in GSM and the constant modulus property. The 
motivation for these approaches is the belief that minimization of a convex linear 
combination of the two cost functions in Equation (5) and Equation (7) will yield 
better performance than minimization of either cost function alone, since the two 
cost functions cannot be simultaneously minimized in general. The resulting cost 
function is given by 
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wfy(0-^*-i(0 ) +(i-Y)(w?y(0-^o(0 



H, 



(9) 



As in the conventionally training adaptive equalizer for use in GSM systems, Tq is 
the time interval occupied by the embedded training signal (FIG. 20). The time 
period-dFfOver which the CM algorithm is applied can be either the entire time 
slot, or the time slot excluding the training signal. Either way, the closed-form 
solution is given by 



w,=[YR^y5^ + (l-Y)R 



(7b) 

yy 



-1 



yRS1,+(1-y)rS^ 



(10) 



The convex linear combiner weight y controls the influence of the two cost 

functions and thus also controls the emphasis given to the two distinct types of 
10 knowledge used in them. Values of y closer to unity give more influence to the 

exploitation of the CM property, whereas values closer to zero give more 

influence to the exploitation of the known training signal. 

It is noted that a special case of the TACM arises when the initial condition 

is zero (i.e., w = 0) and y = 1 . In this case, the TACM yields an algorithm in 
15 which the known training signal is used to initialize the conventional CM 



algorithm at iteration A: = 1 with = Ryy°^] R 



^y^o ' 



In the event that SQ(t) is a constant-modulus signal (e.g., as it is for GMSK, 



CPFSK, and other MSK signals), then SQ(t) should be pre-normalized to have unity 
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modulus prior to its use in the TACM algorithm. In the event that SQ(t) is not a 
constant-modulus signal, then an adaptive normalization constant must be inserted 
into the TACM cost function above to yield the adaptively normalized TACM 
(AN-TACM) algorithm: 



wf y(0 - (0 ' ) + (1 - Y )(|wf y(0 - cs^ {tA . 

ItGTi \ I teTQ 



(11) 



where c is a complex scalar. The solution to this problem is obtained by equating 



y^he complex gradients with respect to and c to obtain 



Y R^^'^ + (1 - Y )Rt?1w, - (1 - Y )RS^c = Y r£> 



y^o 



^y^t-i 



Soy *^^o^o^ 
10 This system of equations is solved by 



(12) 

(13) 







"Rii 


Rl2 


-1 




c 




.^21 


^22. 




0 



15 



(14) 



where Rii = YR^y5^ + (l-Y)R^yJ'^\ R,2= -(I-Y)rS^ R2i = RS^ai^^ 

R22 ~ ~ RiQ^o ■ 

Training-Constrained CM TTCCM) 
Partially Blind Adaptiv e Equalization 
The TCCM differs from the TACM in the way that it attempts to exploit 
both types of knowledge. In particular, when the number of filter coefficients is 
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greater than or equal to the number of samples of the training signal, the linear 
system used by the conventionally trained adaptive equalizer (implied in Equation 
(5)) is under-determined. That is, the training signal is sufficient only to restrict 
Wk to lie in some subspace but not to define it completely. Specifically, in this 
5 case all that can be said so far is that yv^ is given by 



in the null space of Ryy^ ? and defined in the obvious way. The 

TJ 

constraint on ensures that y(/) exactly reproduces the training signal over 

10 time interval Tq. The CM property is then used to select y^, where k denotes as 
before the iteration in the block least-squares implementation of the CM 
algorithm. The resulting cost function is obtained simply by substituting Equation 
(15) into Equation (7) and minimizing with respect to Vj^ to obtain the solution 




(15) 



where superscript ( Y denotes the Moore-Penrose pseudoinverse, Vj^ is any vector 




(16) 



15 where 




(17) 



and the columns of V„ are the eigenvectors associated with the negligible 



49 



EM519050807US 



• 



eigenvalues of Ryy^ • 

_ An alternative algorithm can be obtained by first solving for the CM weight 
vector according to Equation (8) and then projecting this vector onto the affine 

space jwj^ = w^^^^ + V„c^:for all . The resulting solution is 

= W^^O) + V„ Vf (R^yy'^)"' RS^_, (18) 

TJ 

and thus y{t) exactly reproduces the training signal over the time interval Tq. 
As with the TACM, if the desired signal S(,(t) has constant modulus, as is the case 
with OMSK, CPFSK, and other MSK signals, then So(t) should be normalized to 
have unity modulus. Finally, each of the two versions of the TCCM algorithm 
may be extended, through the adaptive normalization technique used in deriving 
the AN-TACM algorithm, to accommodate desired signals having non-constant 
modulus. These extensions are straightforward to obtain and so are not elaborated 
upon here. 

Additional Adaptation Methods 
In addition to the aforementioned TACM and TCCM methods and their 
various extensions and special cases, other adaptation methods can also be used to 
select the values of the coefficients in the filter. For example, the so-called 
decision-direction method of adaptation, in which a first phase of adaptation 
proceeds under control of some method (e.g., any of the aforementioned methods) 



50 



EM519050807US 



• 



to produce a medium-to-high quality signal estimate which can then be 
demodulated, and the resulting detected bits are used to form a known training 
signal which can be used together with the original perfect training data in a 
second phase of adaptation (e.g., through the conventional least-squares method). 
5 E. Preferred Embodiment of an Overall Receiver 

FIG. 22 depicts a preferred embodiment of an overall receiver apparatus 36 
in accordance with the present invention. RF signals at the outputs of the M 
antennas 38 are coherently downconverted to complex baseband and digitized by 
a Coherent M- Antenna Collection Apparatus 40, the output of which is a digitized 
^ 10 discrete-time sequence^ x 1 comple^^^^veeter, denoted elsewhere herein as x(t). 
This vector of received data is then provided as input to a Filter Apparatus 42, 
which takes as its other input the filter coefficients (W) provided by an Adaptation 
Apparatus 44, to produce a digitized discrete-time sequence oidxl complex 

vector signal estimates s^(/7) . During the adaptation process, a switch 46 at the 
15 output of the Filter Apparatus 42 directs the output of the Filter Apparatus 42 into 
the Adaptation Apparatus 44; upon the detection of convergence by the 
Adaptation Apparatus 44, such as when a prescribed number of iterations are 
completed, the switch 46 is moved to direct the output of the Filter Apparatus 42, 
the final signal estimate s{n) , to subsequent processors such as demodulators that 
20 will render bit decisions to recover the binary data streams conveyed by the 
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desired signals. An intermediate output of the Filter Apparatus 42 (denoted by y(t) 
herein) is also provided as input to the Adaptation Apparatus 44. 

A preferred embodiment of the Coherent M- Antenna Collection Apparatus 
40 is shown in FIG. 23. This apparatus is controlled by user-supplied values for 
the RF tuning frequency and the A/D sampling clock frequency to control the 
local oscillators 48 used for the downconversion from real RF to complex 
baseband (also known as in-phase and quadrature-phase signals) and for the A/D 
converters 50, respectively. The RF front-ends 52 that perform this 
downconversion can be implemented in a variety of ways known well to those 
skilled in the art of RF design. Preferably the sampling clock frequency of the 
A/D converters is equal to an integer multiple of the bit rate^J,. With the exception 
of the LCL and RI filters for which a sampling rate offf^ is needed, 2yj, represents a 
good compromise between the temporal resolution provided by oversampling and 
the excess of data and associated computational load induced by oversampling. 
The digitized discrete-time inphase and quadrature signals at the output of each 
pair of A/D converters are combined to form a representation of a single complex- 
valued digitized discrete-time signal. The M signals so formed, one per chain of 
hardware associated with a single antenna, are grouped into a digitized discrete- 
time Mxl complex vector of received data. 

A preferred embodiment of the Filtering Apparatus 42 is shown in FIG. 24, 
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The. Filtering Apparatus 42 accepts as input the Mx 1 complex vector of received 
data and the filter coefficients output from the Adaptation Apparatus 44. The 
received data is first processed by a Frequency-Shifting, Time-Shifting, and 
Polarized Switching (i.e., conjugation, real-part selection or imaginary-part 
selection) Apparatus 54 whose purpose is to form the vector y(t). The record of 
such transformed data y(t) is stored in a Signal Buffer 56, which can be 
implemented using dynamic RAM chips or other high-speed storage media, to 
allow y(t) to be made available to the Adaptation Apparatus 44 over multiple 
iterations of the adaptation algorithm implemented therein. A Matrix- Vector 
Multiplier 58 applies the filter coefficients to the received data, producing output 

s^(«) for use in the subsequent iteration (number k+ I) performed by the 
Adaptation Apparatus 44. 

A preferred embodiment of the Frequency-Shifting, Time-Shifting, and 
Polarized Switching Apparatus 54 is shown in FIG. 25 for the particular filtering 
structure corresponding to the 3 -path LCL-FRESH filter described previously 
herein. It will be appreciated that the other filters previously described herein 
could be alternatively used. The apparatus accepts as an input the Mx 1 vector 
x(t) of received data and produces as an output the 3LMx 1 vQctor y(t) according 
to Equation (4). The Tapped Delay Lines (TDL) 60 produce sets of L time-shifted 
versions of their respective inputs. The Polarized Switcher 62 which, in this 
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embodiment is a conjugator, and the conjugator 62 and Frequency-Shifters 64 are 
depicted using standard block symbols from signal processing and 
communications. Similarly, the preferred embodiments of the Frequency-Shifting, 
Time-Shifting, and Polarized Switching Apparatus 54 for the particular filtering 
structures corresponding to the LCL-FRESH-FSE (2 alternate embodiments), the 
LCL-FSE (2 alternate embodiments), the RI-FSE (a second alternate embodiment 
can be obtained for the RI-FSE, similar to the LCL-FRESH-FSE and LCL-FSE), 
the LCL, and the RI filters are shown in FIG. 26 through FIG. 32, respectively, in 
which the intermediate signal vector at the output is 4LMx 1 or 2LMx 1 . Note 
that the embedded training sequences that appear in the desired signal(s) must also 
be bit-rate sampled by a Bit-Rate Sampler 66 and frequency-shifted by /* prior to 
using them in adaptation methods applied to the LCL-FSE, RI-FSE, LCL, and RI 
filters; it is similarly noted that this frequency shift must be compensated for at 
the filter output by a conjugate frequency shift / if it is desired to obtain an 
unshifted estimate of the desired signal. However, many GMSK demodulators 
require their inputs to be shifted by /* anyway, and so the conjugate frequency 
shift / and the subsequent frequency shift /* cancel each other. It is also noted 
that the RI-FSE and RI filters can be more efficiently implemented due to their use 
of real-valued intermediate signal vectors (and consequent arithmetic 
simplification of filter application and adaptation), instead of the complex-valued 
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intermediate signal vectors that are used in the LCL-FSE and LCL filters. 

Referring to FIG. 33, Polarized Switcher 62 is depicted in a functional block 
diagram form. It can be seen at this point that Polarized Switcher 62 represents a 
family of filter design parameters wherein the real 68 and imaginary 70 parts of 

5 the complex input signal are optionally selected or deselected by selectors 72 
according to the conjugate cycle frequencies and input/output sampling rates 
shown in Table 1 . Note also that the sign of the imaginary part of the input signal, 
if selected, is either unaffected or flipped by switch 74. In other words, the 
complex input signal is either conjugated, its real part is selected, or its imaginary 

10 part is selected. It will be appreciated, however, that Polarized Switcher 62 is 
preferably implemented as a single function conjugator or real-part selector or 
imaginary-part selector rather than an apparatus that provides for the altemative 
selection of each of those functions. 



15 (TDL) 60 is shown. This is a standard implementation of what is known in the art. 
The LMx 1 output vector b(t) produced thereby is defined simply for any Mx 1 
input vector a(t) to be 



Referring to FIG. 34, a preferred embodiment of the Tapped Delay Line 



a(0 
a(f-l) 
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Thus, TDL 60 can accommodate either a real input or a complex input, with 
appropriate modification to the underlying unit-sample delay devices and signal 
paths. 

Referring now to FIG. 35 and FIG. 36, alternative embodiments of 
Adaptation Apparatus 44 are shown where TACM and TCCM are the chosen 
adaptation algorithms, respectively, and all of the d desired signals in the x 1 
vector SQ(t) have constant unit modulus. Adaptation Apparatus 44 takes as inputs 
the intermediate signal datayft), provided by the signal buffer 56 in the Filtering 
Apparatus 42, and the output signal s^(0 computed by the Matrix- Vector 
Multiplier 58 in the Filtering Apparatus 42. It provides as output the matrix of 
filter coefficients, and it also controls the output switch 46 shown in FIG. 22. 
Adaptation Apparatus 44 further moves the switch 46 from its lower position, in 
which Filtering Apparatus 42 output is directed to the Adaptation Apparatus 44, to 
its upper position, in which the Filtering Apparatus 42 output is the <3f x 1 vector 

s(^) of desired signal estimates, suitable for subsequent use by an appropriate 
demodulator. It should be understood by those skilled in the art of implementing 
(e.g., on programmable digital signal processing chips) that alternate embodiments 
of the Adaptation Apparatus 44 follow directly from the descriptions of the 
various versions of TACM and TCCM in this document. 

As can be seen, the present invention can be used to extract a signal of 
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interest from a plurality of spectrally overlapping communications signals, 
separate and remove distortion from interfering co-channel signals and suppress 
adjacent-channel interfering signals of the GMSK or MSK type. These signals 
characteristically have real and imaginary components, exhibit spectral and 

5 temporal overlap, exhibit temporal redundancy, have conjugate cycle frequencies 
equal to twice their carrier frequencies plus and minus one-half of their data bit 
rate, and exhibit conjugate spectral redundancy for spectral components having 
frequencies separated by their conjugate cycle frequencies. The present invention 
provides for time-shifting, frequency-shifting, and conjugation or extraction of 

10 real or imaginary parts of complex signals, to produce an estimate of the signal of 
interest. While the description above contains many specificities, these should not 
be construed as limiting the scope of the invention but as merely providing 
illustrations of some of the presently preferred embodiments of this invention. 
Thus the scope of this invention should be determined by the appended claims and 

15 their legal equivalents. 
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